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ABSTRACT 



I he purpose of this thesis is to find suitable ways to design feedback, compensators 
for liigh order systems by using Root-Locus methods. 

As a starting point we will examine a motor amplidyne system and a position con- 
trol system that were previously designed using Bode methods. I'hen we generalize the 
method and extend it to other systems. 

1 he final subject of this theiis is to design feedback compensators as filters by using 
state feedback coefficients to define zeros of the filter, then we extend this idea to build 
cascade filters. 
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I. INTRODUCTION 



Modern control systems may be compensated by placing a suitable filter in either 
the forward (cascade) path or in a feedback path. 

Mainly feedback systems have the advantage that the output follows the conunand 
more accurately, so the steady state error is less, the bandwidth is greater, and the speed 
of response is faster. Also we can say that the elfcct of external disturbances and internal 
changes in the parameters or structure of the system is significantly reduced. 

The conventional and more conmion way to design a Feedback compensator to 
meet a given set of specifications is the Bode plot method, by determining a suitable gain 
cross-over frequency and 1/H curve. 

In this thesis we will try to find some rules to do the same design by using Root- 
Locus methods. To achieve this, first we will examine several designs then try to put the 
general ideas in a set of rules. 

Simulation studies employed Dynamic Simulation Language (DSL) and Ewald 
package in the IBM 3033/4381 main frame. 
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II. PRELIMINARY EXAMPLES 



A. HIGH GAIN MOTOR AMPLIDYNE SYSTEM 



To study the ideas for feedback compensation first we choose an industrial design 
that exists and is available for comparison. Figure 1 shows the whole compensated 
system, where 




Figure 1. Original Compensated System 



A', = 14.8 


(1) 


. 1209.8 

s(s-f57.96X5 + 0.12+J9.32) 


(2) 


0.137sYs + 62.5)(s-f 71.65) 
“ (s + l)(s + 5.26)(s + 625) 


(3) 



Figure 2 shows the open loop Bode plot for the uncompensated system. The un- 
compensated system is badly unstable. 
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We may use two difTerent kinds of block diagram manipulation to analyze this sys- 
tem. First is the ordinary manipulation shown at Figure 3, where 




Figure 3. Ordinary Manipulation 



/f, = 14.8 (4) 

6096553.14(.r -I- l)(s -I- 5.26)(s -I- 625) 

^ s{5 -I- 0.003+/).066)(j -I- 59.35)(s -I- 79.2)(5 -t- 275.45+J845.8) 

y\nalysis of this transfer function provides the open loop Bode plot, Root- Locus, 
magnified Root-Locus and time response of the system which are given by Figures 4 - 
7. Since there is no pole in the right half plane, we can get phase margin and gain margin 
from these plots. 
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Figure 6. Comp Sjsj teiii Magnified Root-Locus 
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Figure 7. Comp System Time response 
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1 he second type of manipulation is done to get more information using tlie Root- 
Locus method, 'flie first step lor this is shown at Figure 8. Ihen with the appropriate 
block diagram reduction the system becomes as in Figure 9, where 




Figure 8. First Step For Manipulation 




Figure 9. Final System For Root-Locus Examination 



^ 17905 

{s -F 579.97)(s -F 56.14)(s - 26.98+y45.19) 



( 6 ) 
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0.009265^(5 + 62. 1 6)(5 + 71.92) 

“ (5 + 1)(5 + 5 . 24)(5 + 625.33) 

By using the loop transfer function GII(s) to draw the Root-Locus we may examine 
the system for root movements. Since the G(s) has two poles in the right half plane, the 
Bode plot cannot be interpreted from the point of view of gain margin or phase margin, 
but it shows the stability of the system. 

The root locus for the uncompensated system is shown at Figure 10. 

From this point we can proceed to explore the design procedure of the designer. 
First of all, since this design was done using Bode design methods, it will be helpful to 
get the G1 and 1 ,'H magnitude curve which is shown at Figure 1 1 for the original system 
at Figure 1. 

As we may see from the block diagram there are zeros at the origin. For this system 
we need at least two zeros to keep the system type number unchanged. Also the spec- 
ifications show that we want to have a very high error coefficient for steady state accu- 
racy and a very narrow bandwidth. Because of these characteristics, there is a dipole near 
the origin. It is possible to see this on the Bode plot of Figure 1 1. 

The designer could have put the third zero ver}' close to origin but in that case 
building the system might difficult. 

Since other poles are far away compared with these zeros and dipole, these compo- 
nents u ill decide the system behaviour. 

fhe Bode plot, Root-Locus and magnified Root-Locus around the origin for the 
loop transfer function GIl(s) are given by Figures 12 - 14. Since we have roots in the 
right half plane, the Bode plot cannot be interpreted for phase and gain margin values. 

As a result of this problem we may say that for high error coefficient and narrow 
bandwidth systems we may need zeros at the origin and a dipole very close to the origin. 
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Figure 10. Uncompensated System Root-Locus 
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Figure 11. G and 1/H Magnitude Curves 
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Figure 12. Bode Plot For Loop Transfer Function GH(s) 
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Figure 13. Root-Locus For Loop Transfer Function GH(s) 
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Figure 14. Magnified Root-Locus For Loop Transfer Function GH(s) 
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B. POSITION CONTROL SYSTEM 



As a second example we chose a position control system having an amplidyne and 
a DC armature controlled motor as its power element. Figure 15 shows the original 
compensated system, where 




Figure 15. Original Compensated System 



K , = 0.233 


(8) 


A'/, = 0.061 


( 9 ) 


P 12000 

' 0.0785+1 


(10) 


r ’ 

^ 1.885(0.04025+1) 


(11) 


Cl 

II 

p 


(12) 


G4 = 0.00465^ 


(13) 


2 




^5- , . 


(14) 
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Figure 16 siiows the open loop Bode plot for the uncompensated system. I he un- 
compensated system is badly unstable. 
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Figure 16. Uncoiiipensated System Open Loop Bode Plot 
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Again we may use two difTerent kinds of block diagram manipulation to analyze this 



system. 

First is the ordinary manipulation shown at Figure 17, where 




Figure 17. Ordinary Manipulation 



a; = 0.233 

12000(.s^ + 6.5s+ 10) 

0.00595^ + 3.628555** + 247.1445s^ + 14.45s^ + 18.8s 



(16) 



Analysis of this transfer function provides the open loop Bode plot, Root-Locus, 
magnified Root-Locus, and time response of the system which are given by Figures 18 
- 21 

Since there is no root in the right half plane, we can get phase and gain margins from 
this plot. 
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Figure 18 , Compensated System Open Loop Bode Plot 
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Figure 19. Compensated System Root-Locus 
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Figure 21. Compensated System Time Response 
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1 he second type of manipulation is done to get more information using the Root- 
Locus method. The first step for this is shown at Figure 22. 




Figure 22. First Step For Manipulation 



Then with the appropriate block diagram reduction the system becomes as in 
Figure 23, where 




Figure 23. Fihal System For Root-Locus Examination 
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(17) 



473898.3 

(5- 26.72+y66.9j(5 + 91.74) 

^ 0.0()12.v''(:? + 72.39) 

(5 + 2.5)(5 + 4) 



(18) 



By using the loop transfer function GH(s) to draw the Root-Locus we may examine 
the system for root movements. Since G(s) lias two roots in the right half plane, the 
Bode plot of this cascade combination cannot be interpreted from the point of view of 
gain margin or phase margin, but it shows the stability of the system. 

The root locus for the uncompensated system is shown in Figure 24. From this 
point we can proceed to explore the design procedure of the designer. 

First of all, since this design was done using Bode design methods, it will be helpful 
to get the G and 1,'H magnitude curve which are shown in Figure 25 for the original 
system in Figure 15. 

As we may see from the block diagram there are three zeros at the origin. As in the 
first example we have to have at least one zero at the origin to keep the system type 
number unchanged. Also since we want to keep the error coefficient unchanged we have 
to have one additional zero at the origin. 1 he reason for the third zero is to provide a 
dipole near the origin to make the system have a very high error coeflicient for steady 
state accuracy and a very narrow bandwidth. 

For the Gll(sj, the Root-Locus and a close look around the origin to the Root- 
Locus are given in Figures 26 and 27. 

Now wc may want to see the effect of removing some of the components from the 
compensator. First, if we remove the third zero the Root-Locus becomes as in 
Figure 28. The system is unstable. Now we may want to bend the loci by moving one 
of the poles close to the origin toward left. If we start to move the pole which belongs 
to the dipole at s = -2.5; Figure 29 shows the condition when this pole is at s = -50. 
The system is still unstable. By putting this pole further from the zero at s = -72.39, 
we can bend the loci toward the left half plane. We may interpret this as follows : when 
we move the pole to the left we increase the sum of the poles and therefore the sum of 
the roots. This affects the Root-Loci by moving the asymptote centroid to left, so the 
loci from the right half plane cross into the left half plane. Figure 30 shows the Root- 
Loci for this condition. Now we may want to e.xamine the specifications of the system 
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and compare them with our specifications. Since the dominant roots are far away from 
the origin the accuracy and damping will decrease, and the bandwidth of the system will 
increase. These do not agree with our specifications. 

After all these trials and sample designs we may summarize the results as a set of 
rules for design of feedback compensation : 

1. To keep the system type number unchanged put at least the same number of zeros 
at the origin as the original system's type number. 

2. lo keep both the system type number and the error coellicient same unchanged, 
the number of zeros at the origin should be N + I, where N is the type number 
of the uncompensated system. 

3. fo have high error coeflicient and narrow bandwidth for steady state accuracy in- 
clude a dipole close to the origin. Put the zero of the dipole at the origin for ease 
in physical realization. 

In the next chapter we will apply these rules to several different kinds of problems 
and see whether they work or not. 



25 



iMflG RXIS 

-100. C -7S.0 -50.0 -25.0 0.0 25.0 50. C 75.0 lOO.C 















« 

• 

/ 


m 

0 














/ 
























' 












7\ 












































\ 
















\ 

$ 


» 

0 



t 



-100.0 -BO.O -60.0 -40.0 -20.0 0.0 20.0 40.0 CO.) 



REAL AXIS 



Figure 24. Root-Locus For Uncompensated System 




Figure 25. G and 1/H Magnitude Curves 
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Figure 26. Root-Locus For The Loop Transfer Function GH(s) 
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Figure 28. Tlie S.inie System With Tuo Zeros .it Origin. 
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III. DESIGNING COMPENSATORS 



A. NON-MINIMUM PHASE SYSTEMS 

Before studying systems that have all poles and zeros in the left half plane, we may- 
want to examine non-ntinimum phase systems. If we have a root or roots in the right 
half plane we are not able to use Bode methods for designing compensators, therefore 
the Root-Locus method will be the only useful tool. 

To illustrate we chose the following transfer function, which is the roll mode of a 
vertical take-olf aircraft. Figure 31 shows the original system, where 




. 77.187(5 -f 0.225 +j0.6607) 

(5 - 0.0434j(5 -I- 1.464X5 4- 0.21 +y0.844j(5 -I- 48.35 ±J49.34j ^ ’ 

As we may see from the transfer function, there is a real pole in the right half plane. 
For this reason we can not use Bode methods to compensate the system. The open loop 
Bode plot, Root-Locus and magnified Root-Locus for this system are given by Figures 
32 - 34. 

To find the closed loop roots we make the block diagram manipulations. Then the 
closed loop transfer function becomes : 



77.187(5 -h 0.225 +y0.6607) 

“ (5 - 0.036)(5 4- 0.2 1 + j 0.84)(5 4- 1 .474X5 4- 48.34 ±y48.72) 



( 20 ) 
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Figure 32. Uncompensated System Open Loop Bode 
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Figure 33. Root-Locus For Uncompensated System 
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From this point on we may start to think in terms of Root-Locus and loop transfer 
functions. 

Since we will want to have loci from the pole in the right half plane toward the left, 
wc have to put cither a zero or a pole to the left of this pole. Putting a zero at the origin 
will not help us, because in that case we always will have a root in the right half plane. 
So we should put the compensator poles and zeros somewhere in the left half plane. 
Choosing the specific zero and pole locations is an engineering decision and mostly de- 
pends on the specifications. With this problem we will try to illustrate three dilferent 
combinations. 



First put a zero close to the origin and a pole far away, both in the left half plane. 
In this case the loop transfer function becomes as in Figure 35, where 




G,(s) 



77.187(5 + 0.225 T;0.6607) 

(5 - 0.036)(s + 0.21 + y0.84)(5 + 1.474Xs + 48.34 + J48.72) 



( 21 ) 



//,( 5 ) 



+ 1 ) 

(5 + 50) 



( 22 ) 



The Root-Locus and magnified Root-Locus for the loop transfer function are given 
in Figures 36 and 37. As we may see from the plots there is a root-locus segment be- 
tween the pole in the right half plane and the zero at -1. 

To find the value for A',, we should get the tabulated points for the Root-Locus and 
see for which gain all of the roots are in the left half plane. For this particular filter the 
minimum value for A', was 203.7. We selected A', = 500 arbitrarily. For this value, our 
system becomes as in Figure 38, where 
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Figure 36. Root-Locus For Loop Transfer Function 
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Figure 37. Magnified Root-Locus For Loop Transfer Function 



77.187(5 + 0.225 T7'0.6607) 

“ (s - 0.0434)(5 + 1.464)(5 + 0.21 + j0.844j(5 + 48.35 + J49.34) 



(23) 



//,( 5 ) 



500(5+ 1) 
(5 + 50) 



(24) 




Figure 38, Compensated System 



If we make the block diagram manipulation, we get the system shown in 
Figure 39, where 




Figure 39. Resulting Equivalent G 

77.187(5 + 0.225+j0.6617)(5 + 50) 

(5 + 0.0269)(5 + 0.235+y0.842)(5 + 1.535)(5 + 49.67)(5 +48.41+y48.63) 

The Root-Locus, magnified Root-Locus, open loop Bode plot and time response for 
the compensated system are given in Figures 40 - 43. As we may see from the Bode plot 
the system is stable but we do not have enough gain. 
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Figure 41. M.ignifed Root-Locus For Final System 
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Figure 42. Bode Plot Without Outside Gain 
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Figure 43. Time Response Without Outside Gain 
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We may handle this by putting another gain block outside the minor loop. I he final 
system then becomes as in Figure 44. 




Figure 44. Final Compensated System 



K may be selected according to the specifications. We chose K = 100 for illus- 
tration. 






77.187(5-f 0.225 +;0.6607) 

{s - 0.0434)(5 + 1.404)(5 + 0.21 + J0.844)(5 + 48.35 + J49.34) 



(26) 



//,( 5 ) 



500(5-1- 1) 
(5 -f- 50) 



(27) 



The open loop Bode plot and time response for this particular values are given in 
Figures 45 and 46. We may change the loop gain according to the specifications. Flfects 
of changing the loop gain will be shown for the third combination. 
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Figure 45. Bode Plot For Final System Loop Gain = 100 



46 



^ o 

rvi 




3 



Figure 46. Final System Loop Gain = 100 
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The second type of configuration puts two complex zeros close to the origin and two 
poles far away, all in the left half plane. In this case the loop transfer function becomes 
as in Figure 47, where 
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Figure 47. Loop Transfer Function Form 



G,(5) = 



77.187(5 + 0.225 +70.6607) 



(5 - 0.036)(5 + 0.21 + 70.84)(5 + 1.474J(5 + 48.34 + 748.72) 



7/2(5) = 



^2(5 + 1 +7) 

(5 + 20)(5 + 30) 



(28) 

(29) 



The Root-Locus and magnified Root-Locus for the loop transfer function are given 
in Figures 48 and 49. As we can see from the plots we moved the loci toward the left 
half plane. 
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Figure 49. Magnified Root-Locus For Loop Transfer Function 
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To find a value for , we should get the tabulated points for the Root-Locus and 
see for which gain all of the roots are in the left half plane. For this particular filter the 
iiiinimuin value for Aj was 2125.8. We selected K2 = 2500 arbitrarily. For this value, our 
system becomes as in Figure 50, where 

. . ^ 77.187(5 -H 0.225 +;0.6607) 

(s - 0.0434)(s + 1.464)(s -t- 0.21 + ;0.844)(s -t- 48.35 + ^49.34) 



. _ 2500(5 + 1 +j) 
“ (5 + 20)(5 + 30) 



(31) 




Figure 50, Compensated System 



If we make the block diagram manipulation, we get the system shown in 
Figure 51, where 




Figure 51. Resulting Equivalent G 



77. 1 87(5 + 0.225 + jO.661 7)(5 + 20)(5 + 30) 

“ (5 + 0.007)(5 + 0.22+;0.84)(5 + 1.427)(5 -f- 24.5+j6.3)(5 -f- 48.87+y48.3) 
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The Root-Locus, magnified Root-Locus, open loop Bode plot and time response for 
the compensated system are given in Figures 52 - 55. As we may see from the Bode plot 
the system is stable but the gain is too low. 
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Figure 53. Magnified Root-Locus For Final System 
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Figure 54. Bode Plot Without Outside Gain 




Figure 55. Time Response Without Outside Gain 
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We may handle this by putting another gain block outside the minor loop. 1 he linal 
system then becomes as in Figure 56. 




Figure 56. Final Compensated System 



K may be selected according to the specifications. We chose K = 100 for illus- 
tration. 



77.187(5-1-0.225 TjO.6607) 

~ [s- 0 . 0434)(5 - t - 1 . 464)(5 - f - 0.21 + y 0 . 844)(5 4 - 48.35 + J 49 . 34 ) 



(33) 



2500 ( 5 + 1 +;) 
(5 + 20)(5 - t - 30 ) 



( 34 ) 



1 he open loop Bode plot and time response for these particular values are given in 
Figures 57 and 58. We may change the loop gain according to the specifications. Efiects 
of changing the loop gain will be shown with the next combination. 

The third type of combination puts two complex zeros and one real pole close to the 
origin and one pole far away, all in the left half plane. In this case the loop transfer 
function becomes as in Figure 59, where 
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Figure 57. Bode Plot For FinnI System Loop Gain = 100 
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Figure 59. Loop Transfer Function Form 



77.187(5 + 0.225 TyO.6607) 

^ (s - 0.036)(5 + 0.21 + yo.84)(5 + 1.474)(5 + 48.34 + J48.72) 



(35) 



A3(5 + 2+y0.5) 
(5+1)(5 + 50) 



(36) 



1 he Root-Locus and magnified Root-Locus for the loop transfer function are given 
in Figures 60 and 61. As we may see from the plots the locus is moved toward the left 
half plane. 
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Figure 60. Root-Locus For Loop Transfer Function 
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To find a value for Aj, we should get the tabulated points for the Root-Locus and 
see for which gain all of the roots are in the left half plane. For this particular filter the 
niiniinum value for K\ was 81.6. We selected A'j = 250 arbitrarily. By this, our system 
becomes as in Figure 62, where 




Figure 62, Compensated System 



77.187(5 4- 0.225 T/).6607) 

(5 - 0.0434)(5 4- 1.464)(5 4-0.21 T;0.844)(5 4-48.35 +y49.34) 



(37) 



250(5 + 2 4- 0.5J) 
( 5 + 1)(5 4- 50) 



(38) 



If we make the block diagram manipulation, we get the system shown in 
Figure 63, where 






77.187(5 4- 0.225 + J0.6617)(5 4- 1)(5 4- 50) 

(5 4- 0.19)(5 4- 0.236 + y0.87)(5 4- 1.52)(5 -f 49.84)(5 4- 48.37 + 7‘48.6) 



( 39 ) 



c 

► 




Figure 63. Resulting Equivalent G 



The open loop Bode plot and time response for the compensated system are given 
in Figures 64 and 65. As we can see from the Bode plot the system is stable but the gain 
is too low. 
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Figure 64. Bode Plot Without Outside Gain 
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Figure 65. Time Response Without Outside Gain 
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We may handle this by putting another gain block outside the minor loop. The final 
system then becomes as in Figure 66. 




Figure 66. Final Compensated System 

K may selected according to the specifications. We chose difierent values for illus- 
tration. 



G{s) = 



77.187(5-1- 0.225 +;0.6607) 

(s - 0.0434)(s -b 1.464)(5 4- 0.21 T J0.844){5 + 48.35 + ^'49.34) 



(40) 



(5 -b 1X5 + 50) 



(41) 



The open loop Bode plot and time response for K= 100 are given in Figures 67 and 

68 . 
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Figure 67. Bode Plot For Final System Loop Gain = 100 



67 



Olf'i 




"5 



Figure 68. Final System Loop Gain = 100 
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We can change the loop gain according to the specifications. Ihgures 69 and 70 show 
the time responses for this system when this gain is 500 and 1000. As we may see from 
the time response we can change the amount of overshoot by changing this gain. 

As a result of this problem we can sunuuarize the steps for designing feedback 
compensators for non-minimum phase systems as following : 

1. Find the closed loop transfer function and get the Root-Locus for this. 

2. Pick the filter poles and zeros according to the Root-Locus rules. 

3. Get the Root-Locus for the loop transfer function. Tabulate the gain values, and 
select the gain that puts all the roots to the left half plane. 

4. Rearrange the system. Look for the Bode plot and time response. If they are ac- 
ceptable, leave it that way otherwise select another gain outside the minor loop. 
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Figure 69. Final System Loop Gain = 500 
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Figure 70. Final System Loop Gain = 1000 
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B. FLEXIBLE ROCKET CONTROL SYSTEM 



As a second example we chose the control system of a flexible rocket which is shown 
in Figure 71. This example was extracted from Ref. 1. 



Structure Structure and 

Demod. Filter Servo Rigid Body 




Figure 71. Original Uncompensated System 



is the transfer function for the rigid body and is defined as 



Gr{s) = 



7.21 



[s+ 1.6)(5- 1.48) 






333 



(5 + 333) 



(42) 

(43) 



^ . 2750 



5^ + 42.25 + 2750 



(44) 



(5^ + 705 + 4000)(5^ + 225+ 1 2800) 

^5^7(5) — ^ , 

( 5 ^ + 305 + 5810)(5^ + 305 + 12800) 



(45) 



And the structure is defined as ; 



■j s 0.686(5 + 53)(5 - 53)(5^ - 152.2s + 14500)(5^ + 153.85 + 14500) 

~ ( 5 ^ + 5 + 605)(5^ + 45.5s + 2660)(5^ + 2.51s + 3900)(s^ + 3.99s + 22980) ^ 
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After performing all calculations the total system becomes as in Figure 72. Since 
the order of the numerator and denominator are quite large, the transfer function is split 
into three parts to fit on the page. 




Go(5) = G„| X G,i X where ; 

39134576.25(5 + 35 + j52.1)(s + 1 1 Tyl 12.6) 

“ (s+ 15+y74.7)(s+ 15+yil2.1)(s + 333)(5 + 21.1 +y4S) 

(s - 21.55 +./21.79)(5 + 72.76 +yi00.24)(5 - 83.1 +y61. 82) 

‘’2“ (5 + 22.75 +y46.29)(s+ 1.255 +y62.44)(s + 1.995 +y 151. 57)(s - 1.48) 



(5+11.6+y24.02) 

(s+ 1.6)(.s + 0.5+y24.6) ^ 

As we may see from the transfer function we have four zeros and a real pole in the 
right half plane. In addition to this most of the poles and zeros which are in the left half 
plane are very close to the imaginary a.xis. 

fhe Root-Locus for the open loop uncompensated system is given by Figure 73. 
The system is unstable. Since there is a real pole in the right half plane we can not easily 
interpret the Bode plot for this system. 

The next step is to find the loop transfer function as in previous examples. By per- 
forming the regular block diagram manipulation the closed loop transfer function be- 
comes G,( 5 ) = G,,|(5) X 6 ^ 2 ( 5 ) X ^,, 3 ( 5 ) which is 



73 




Figure 73. Root-Locus For Uncompensated System 
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39134576.25(5 + 35 +./52.7)(5+ 11 +yll2.6)(5+ 11.6+j24) 
(5+ 333^(5+ 1.9 Tyl51.5)(5+ 15 +jll2.2j(5+ 13.8 +>74.5) 

(5 + 72.8 + >100.2)(5 - 83 + >61.8)(.s - 21.6 + >21.8) 

~ (5 + 1.5 + >60.4)(5 + 19.9 + >52.6)(5 + 1 + >24)(5 - 0.5 + >4.5) 



Gciis) = 



1 

(5 + 24.8 +>44.17) 



( 52 ) 



From this point on we may start to think in terms of Root-Locus design. As we saw 
from the loop transfer function we have a pair of complex poles and four zeros in the 
right half plane. In addition to this most of the poles in the left half plane are very close 
to imaginary axis and there are four zeros in the right half plane. So probably some of 
the roots are going to end up at these zeros. Tliis will give us a gain constraint. 

Since we have poles close to the origin and imaginary axis in the left half plane, 
putting a zero or zeros at the origin would not help us, because in this case we always 
will have roots in the right half plane. So we decided to use a general type of 
compensator which is two zeros close to the origin and two poles far away. The purpose 
of the zeros was to have loci from the poles in the right half plane terminate on these 
zeros. We made several trials with this type configuration but none of them worked. 
An illustration of this is shown in Figure 74, where GXs) is the closed loop transfer 
function and 




Figure 74. Illustration of General Configuration 



98.6(5 + 0.46 + >4.482) 
(5+ 100)(5 + 200) 



( 53 ) 
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The Root-Locus for this case is given in Figure 75, 
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Figure 75. Root-Locus for General Type Compensated System 



As we may see the system is unstable, because roots in tlie left half plane were 
moving faster than the roots in the right half plane. 

With this type of problem, since the system has a 12"' order numerator and 17"' order 
denominator, designing the compensator by just examining the system would be impos- 
sible or require many trials. So we decided to use the function minimization subroutine 
in the DSL package in the mainframe. The parameters used in the subroutine and the 
function minimization program are given in Appendix A and B. 

As the result of the program we found that we have to put the poles of the 
compensator close to the origin and to the left of the zeros which belong to the 
compensator. By putting these zeros and poles in this combination, we make the roots 
in the right half plane move faster. Figure 76 shows the block diagram for the compen- 
sated system, where G,.(s) is the closed loop transfer function and 




Figure 76. First Resulting Compensator of Function Minimization 



7/2(5) = 



14.72(5 -f 0.875 TJ2.058) 
(5 4-2+ y8.343) 



(54) 



1 he Root-Locus for this compensation is given in Figure 77. By getting the tabu- 
lated root locations, we see that for a Root-Locus gain of 0.068 we get all roots in the 
left half plane. Then our final compensator becomes ; 



7/2(5) = 



5 -+ 1.755 + 5 
5^ + 45 + 73.6 



(55) 
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The time response for this case is given in Figure 78. As we may see the system is 
stable but the settling time is a little long. According to the specs we may want a faster 
system than this. 




Figure 77. Root- Locus for First FM Result 
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Figure 78. Time Response for First FM Result 
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Running the Function Minimization program but changing C and vv„ to get an out- 
put with the desired settling time, we got the coeHicients for tliis compensator as shown 
in Figure 79. G,.(s) is the closed loop transfer function and 




Figure 79. Second compensator Resulting from Function Minimization. 



= 



8.11(5+ 1.1531)(5 + 4.878) 
(5 + 3.672 + 75.67) 



(56) 



The Root-Locus for this is given in Figure 80. Again by getting tabulated root lo- 
cations we see that for a Root- Locus gain = 0.124 we get all roots in the left half plane. 
Then our final compensator becomes ; 



7/3(5) = 



5^ 6.03115 + 5.625 
5^ 7.3445 + 45.632 



(57) 



The time response for this case is shown in Figure 81. As we see the settling time 
and amount of overshoot are within specilic limits. 
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Figure 80. Root-Locus for Second FM Result 
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Figure 81 . Time Response for Second FM Result 
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Then our total compensated system becomes as in Figure 82, where G(s) is the Open 
loop transfer function and 




(5+1.1531)(5 + 4.878) 
(S + 3.672 + J5.67) 



(58) 



As a result of this problem we may generalize that, if we have components close to 
the imaginary axis in both left and right half plane, we have to use pole-zero couples 
close to each other as a compensator. In the next chapter building filters by using pole 
placement and state feedback methods will be discussed. 
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IV. DESIGNING FEEDBACK COMPENSATORS USING POLE 

PLACEMENT 



To compensate the systems by state feedback, we may use several computer pro- 
grams to find the feedback coefficients. Sometimes observing or feeding back some of 
the states may not be possible. For such cases we might have to build an estimator to 
estimate the states. 

The purpose of this chapter is to investigate whether it is possible to build a filter 
by using these feedback coefficients as zeros and adding extra poles to the feedback 
path. 

To illustrate this we chose the following plants. 

A. ALL COMPONENTS ARE IN THE LEFT HALF PLANE 
For this case the open loop transfer function of the plant is ; 



G(5) = 



1000 

5(5 -I- 1)(5 + 2)(5 -I- 10) 



(59) 



We want to design a feedback compensator (a filter) to stabilize this system after 
determining the state feedback coefficients. Since the uncompensated system is fourth 
order, in order to have full state feedback we have to name four roots. 

In this thesis we used the SYS (State Variable Systems) (Ref 2) program. 

To use the computer program we have to determine n-1 roots (where n is the system 
order) and find the fourth root from the characteristic equation to enter the program. 
Our desired roots are 



5|_2 — — 1 +y2, 53 — —4 



(60) 



Then H(s) becomes 



H{s) = 5^ + 65^ + 135-1-20 



If we get the GH(s) to find the characteristic equation 



G//(5) = 



1000(5^ -b 65^ -b 135 -b 20) 
5“^ -b 135^ -b 325^ -b 2O5 



(61) 



(62) 
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rile characteristic equation then becomes 

5 " + 10135^ + 60325^ + 130205 + 20000 (63) 

And the roots of the characteristic equation are ; 

5, = -3.988, 52,3 = -0.997 T jl .996, 54 = - 1007.02 (64) 

The root locations obtained with this design are within 0.3 % of the desired values 
and the 4"' root is at —1007.02 . 

Now we may go to the SVS program to find the state feedback coelTicients which 
would place the roots exactly at —4, — 1 T j2, —1007.02. The state feedback coefficients 
we got from SVS are ; 

7 /( 5 ) = 0.99985^ + 5.9825^ + 12.915 + 19.8 (65) 

These coefficients are negligibly difi'erent from the ones we obtained by multiplying 
the desired roots. For accuracy we used the coefficients that we got from SVS. 

For all practical purposes we may accept the which is coelficient of 5 ’, as 1. Our 
state feedback compensated system then becomes as in Figure 83, where G(s) is the 
open loop transfer function and 




H{s) = 5 ^ + 5.9825^ + 12.915 + 19.8 (66) 

H(s) in factored form is 

H{s) = (5 + 0.996 T yi .993)(5 + 3.997) (67) 

Root-Loci for the loop transfer function GIl(s) are given in Figure 84. 
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To preserve the unity feedback we use block diagram manipulation and the system 
then becomes as in Figure 85. 




Figure 85. The Form of The Syslem With Unity Feedback. 



In Figure 85, G(s) is the open loop transfer function and 

//,(5) = s(s + 2.99+ 1.993) (68) 

To convert the feedback compensator to make it more realizable as a filter, we have 
to choose the locations of the extra poles quite far away so they can not alfect the 
transient behaviour of the system. If we put the extra poles into the system form shown 
in Figure 83 we get the required gain easily. The new system then becomes as in 
Figure 86, Where G(s) is the open loop transfer function and 

/f(5 + 3.99)(^ + 0.996 T;1 .993) 

(j + 3U0)(5 + 400)(i + 500) ^ ' 
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Figure 86. System With Extra Poles. 



To find the value of the K we get the tabulated root locations of the Root-Locus for 
the loop transfer function 67/2(5), which is shown in Figure 87. 
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Figure 87. Root-Locus For Filter Feedback Loop Transfer Function. 
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Root-locus gain K to bring the roots where we want them was K= 1.905459 x 10’ 
. By using the same manipulation as in Figure 85 if we redraw the block diagram of the 
system to preserve the unity feedback we get the final compensated system shown in 
Figure 88. 




Figure 88. Final Compensated System 



G(s) is the open loop transfer function and 



G,(s) = 



19.8 X 1.905459 x 10^ 

(5 -t- 300)(5 + 400J(5 + 500) 



(70) 



im = 



1.905459 X 10'.s(5 -t- 2.99 ±jl.99) 
(5 -t- 300)(s -t- 400)(s -t- 500) 



(71) 



The time responses for only state feedback and for the filter compensated case are 
Figures 89 and 90. As we may see from the transient responses since all roots are very 
close to the desired places the time responses are almost e.xactly same. The only disad- 
vantage of this configuration may be to build two identical sets of poles to use in dif- 
ferent places. 

B. NON-MINIMUM PHASE SYSTEMS 

As a second example we chose a non-minimum phase system to check if the idea 
works for both cases. 

For this case the open loop transfer function of the plant is 



G{s) = 



1000 



s(5 - 1)(5 -t- 2)(5 + 20) 



(72) 
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Figure 89. Time Response for Only State Feedback Compensated System 
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Figure 90. Time Response for Filter Compensated S)stem 
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We want to design a feedback compensator (as a filter) to stabilize this system after 
deterntining the state feedback coefficients. Since the uncompensated system is fourth 
order, in order to have full state feedback we have to name four roots. 

To use the computer program we have to determine n-1 roots (where n is the system 
order) and find the fourth root from the characteristic equation to enter the program. 
Our desired roots are ; 



5] 2 — ~1 +Jlt ^3 ~ 

Then H(s) becomes 

ll{s) = 5^ + + 85 + 6 

If we get the GH(s) to find the characteristic equation 

^ 1000(5^ + 55' + 85 + 6) 

GHis) = — ■: ; 

5^ + 215'^ + 185^ — 405 

The characteristic equation then becomes 

5" + 10215^ + 50185^ + 79605 + 6000 

And the roots of the characteristic equation are 

S, = -2.967, 52.3 = -0.988 + yi.Ol, 54 = -1016 



(73) 

(74) 

(75) 

(76) 

(77) 



The root locations obtained with this design are within 1.1 % of the desired values 
and the 4'* root is at -1016. 

Now we may go to the SVS program to find the state feedback coellicients which 
would place the roots exactly at —1 +J1, — 3, — 1016. 

The state feedback coefficients we got from SVS are 

H{s) = 0.99995^ + 4.9955^ + 7.9525 + 5.956 (78) 

These coefficients again are negligibly different from the M(s) we first accepted but 
we used the ones we found from SVS. 

For all practical purposes we may accept the A3 which is coefficient of 5’ as 1. Our 
state feedback compensated system then becomes as in Figure 91. 
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G(s) is the open loop transfer function and 

H(s) = 5^ + 4.995s^ + 7.952s + 5.956 (79) 

Il(s) in factored form is 

//(s) = (s + 0.995 + y0.998)(s + 2.998) (80) 

Ihe Root-Locus for the loop transfer function GM(s) is given in Figure 92. 
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Figure 92. Root-Locus For State Feedback Loop Transfer Function. 
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To preserve unity feedback we may use block diagram manipulation. The next form 
of the system then becomes as in Figure 93 where G(s) is the open loop transfer func- 
tion and 



/y,(5) = 5(5-f 2.494 T 1.316) (81) 




Figure 93. The Form of The System With Unity Feedback. 



To convert the feedback compensator to make it more realizable as a filter, we have 
to choose the locations of the extra poles quite far away so they can not affect the 
transient behaviour of the system. If we put the extra poles into the system form shown 
in Figure 91 we get the required gain easily. The new system then becomes as in 
Figure 94. 
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G(s) is the open loop transfer function and 



K{s + 2.99S)(s + 0.995 + /).998) 
(5 "f 3U0)(5 "f- 400)(5 "f 5U0) 



To find the value of the K we get the tabulated root locations of the Root- Locus for 
the loop transfer function G I { 2 ( 5 ), which is shown in Figure 95. 
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Figure 95. Root-Locus For Filter Feedback Loop Transfer Function. 
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The Root-locus gain K to bring the roots where we want them was 
A' = 1.12571 19 X 10’ . By using the same manipulation as in Figure 93 if we redraw the 
block diagram of the system to preserve the unity feedback we get the final compensated 
system shown in Figure 96. 




Figure 96. Final Compensated System 



G(s) is the open loop transfer function and 

5.956 X 1.12571 19 X lo’ 



G,(s) = 



{s + 300)(s + 400)(5 + 500) 






1.1257119 X 10^5-1- 2.494 ±71.316) 
(5 -f- 300)(s -1- 400)(5 + 500) 



(83) 



(84) 



1 he time responses for state feedback and for the filter compensated case are given 
Figures 97 and 98. As we may see from the transient responses since all roots are very 
close to the desired places the time responses are almost exactly same. The only disad- 
vantage of this configuration may be to build two identical sets of poles to use in dif- 
ferent places. 

C. EXTENDING THE IDEA TO CASCADE COMPENSATION 

The procedure for root placement with the feedback filter used only the loop transfer 
function. We may compensate the system with the same components and same gain by- 
using cascade compensation. If we investigate the systems, the characteristic equations 
are the same with the filter in either the feedback path or in the forward path. Except 
this time instead of the poles, the zeros of the filter are the zeros of the closed loop 
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Figure 97. Time Response For Only State Feedback Compensated System 
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Figure 98. Time Response For Filter Compensated System 
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transfer function but the roots of the closed loop transfer function are in the same places 
as with the feedback compensation scheme. 

To show the result of this idea we chose the same plants. For the first one if we put 
the same filter in the forward path the S3 stem becomes as in Figure 99, where 




Figure 99. Cascade Compensated System 



G(s) = 



icon 

+ I)(s + 2)(s +10) 



H,{s) = 



1.905459 X 10^(s + 1 + J2){s + 4) 
(s + 300)(s + 400)(s + 500) 



(85) 

( 86 ) 



The time response for this system is given in Figure 100. As w-e may see the time 
responses and oscillatory frequencies are almost identical with the feedback compen- 
sation time response which is shown in Figure 90. In this configuration because of the 
derivatives due to the cascade zeros we have faster rise time. 

We checked the roots of the systems for both the feedback and the cascade config- 
uration. I hey were identical within 0.37 *3'o difference. 
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Figure 100. Time Response for Cascade Compensator 
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I he second plant becomes as in Figure 101 if we pul the compensator in tiie for- 
ward patli, where 




Figure 101. Cascade Compensated Second System 



, 1000 

5(5- l)(s-f 2)(5+20) 



7/2(5) = 



1.1257119 X 10 ^( 5 + 1 + j)(5-t- 3) 
( 5 -1-300X5 -f 4U0)(5 -f 500) 



(87) 



( 88 ) 



1 he time response for this case is given in Figure 102. Again as we may see the time 
response is essentially the same as that of the feedback configuration which is shown in 
Figure 98. Here again we have a faster rise time because of derivatives in the cascade 
filter. 

As we may see from the Figures 99 and 100 we put the filter before the plant in both 
cases and we used the DSL (Dynamic Simulation Language) which behaves the same 
way with the hardware. 

Then we decided to simulate the systems with ALCON (Simulation program for 
persona! computers). This program uses arithmethic polynomial calculations to calculate 
the step responses for the systems. 
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Figure 102. Time Response For Cascade Compensator (Second example) 
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1 he resulting time responses from this program for both problems are given in Tig- 
ures 103 and 104. As we see high frequency ripples exist for a short time period which 
we didn't see in the DSL simulation. 

To investigate the reason for this dilference we interchanged the blocks in Figures 
99 and 100, that is, we fed the error signal into the G block and the output from the G 
block into the H block. For this combination resulting step responses are given in Fig- 
ures 105 and 106 which are identical with the time responses we obtained from ALCON. 

For these DSL simulation cases, in the first case the filter is before plant, the plant 
works as Low-Pass filter and the high frequency ripples are not present, but in the sec- 
ond case the filter follows plant, and the ripples appear in the output. 

We could not decide which one of the simulations is true. According to linear the- 
ory, changing the places of the plants should not affect to the behaviour of the system. 

To investigate the reason for this difference is left as a recommendation. 

As a result of these examples we may generalize the design steps as follows : 

1. Select N - 1 desired roots. 

2. Multiply them and find the H(s). 

3. Find GH(s) and get the characteristic equation. 

4. Find the roots of the characteristic equation. If they are acceptably close to the 
desired roots, use H(s) as a state feedback directly, otherwise go to SVS or another 
computer program to find the state feedback coefficients to put the roots exactly 
at the desired places. 

5. Choose the extra poles as far away as possible and put them in the denominator 
ofll(s). 

6. Draw the Root-Locus for the loop transfer function GH(s) (including the extra 
poles) and get the tabulated data. 

7. Select the appropriate gain to have the system roots in the desired places (generally 
this will be very high gain). 

8. If the purpose is to design a cascade filter, put the resulting filter directly in the 
forward path. Otherwise perform the block diagram manipulation as shown in the 
examples and preserve the unity feedback. 
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Figure 103. ALCON Simulation of First Problem 
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Figure 104. ALCON Simulation of Second Problem. 
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Figure 105. Second DSL Simulation of First Problem 
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Figure 106. Second DSL Simulation of Second Problem 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

111 this thesis, the development of a procedure to design the feedback compensators 
by using Root-Locus methods and pole placement methods has been presented. Simu- 
lation results and worked out examples have shown that except for extremely high order 
systems designing the compensators with the Root- Locus method can be used for any 
plant. The most important result is to design them using pole placement concepts and 
to be able to use them either in cascade or in the forward path. 

Simulation results have also shown that the compensator can be determined by a 
Function Minimization subroutine directly from the transfer function of an ideal plant 
based on given specifications. 

B. RECOMMENDATIONS 

1. Further research should investigate methods to satisfy certain specifications while 
designing compensators with the Root-Locus method. 

2. In this thesis, it has been shown that Function Minimization can be used to design 
compensators, but the user should write his own program according to the param- 
eters. A package program may be written for the main frame to make the procedure 
interactive and faster. 

3. Further research may investigate the effect of gain variations on the root locations 
while using a pole placement method so one can use directly multiplication of de- 
sired roots as H(s). 



I 



111 



APPENDIX A. CONSTRAINT PARAMETERS OF FUNCTION 

MINIMIZATION 



Subroutine HOOKE uses the Hooke and Jeeves "direct search" method to find the 
local maximum or minimum of a multi-parameter criterion function, CF. [Ref. 3]. 

1 he algorithm evaluates CF at a base point, X= (^'(1), ..., ^'(A')), then perturbs each 
parameter in turn, by the amount ± STEP{F), and evaluates CF at each new point. 

If none of these points produces a better value of CF, then the stepsizes are de- 
creased by a factor of BETA, and the process is repeated. 

To use HOOKE, the user must initialize the following arguments and, in the main 
program, CALL lIOOKE(X,STEP,N,ITMAX,CFTOL,ALPHA,BETA,CF,Q,QQ, 
W,1PR1NT,M INMAX) 

All of these arguments must be initialized in MAIN, except for X, CF, Q, QQ, and 
W. Recommended values for ALPHA and BETA are ALPHA = 2., BETA = 0.5. All 
of the arrays, i.e., X, STEP, Q, QQ and W, must be declared and dimensioned in .MAIN. 

1 he arguments and their meanings may be explained as following. 

X = The array of N parameter values. The user must supply the initial guesses, ei- 
ther in the DSL program or in MAIN. 

STEP = An array of dimension N containing the initial step sizes to be used in the 
search. 

N = The number of parameters (a positive integer, at most 15). 

ITMAX = The maximum number of function calls to be performed. 

CFTOL = The error in the criterion function to be reached before the program 
terminates ( diflerence between the current value and the previous stage value). 

ALPHA = The factor of (Y - X) which is added to Y to get XNEW; a number 
greater than or equal to 1. 

BETA = The stepsize reduction factor; a number between 0 and 1. 

CF = The value of the criterion function. 

Q, QQ, W = Arrays of dimension N, to be used as work space. They must be de- 
clared and dimensioned in the main program. 

I PRINT = An integer flag. 0 for no intermediate printout, 1 for intermediate 
printout ofX, CF, the number of function evaluation and notification of step-reduction. 
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APPENDIX B, COMPUTER PROGRAM FOR FUNCTION 

MINIMIZATION 



*This program finds the necessary parameters to fit the function to some 
^special transient response. By changing zeta and omgn the standart time 
*response may be generated. 

D COMMON /HAND J/FLAG , ERRFN , A 1 , A2 , B 1 , B2 
TITLE Rocket parameters 

* FUNCTION MINIMIZATION 

*Put initial conditions according to the previous assuming parameters. 
INCON A10=2. ,A20=5. ,B10=1. 4B20=100. 

^Define your arrays to simulate your system 

ARRAY P(1),Q(2),A(1),B(2),D(5),E(5),F(1),G(3),H(9),L(11),K(2),... 

M(2),Y(3),Z(3),S(2),T(4) 

*Put the values for the array coefficients 

TABLE A(l)=333. ,B(1-2)=1. ,333. ,D(1-5)=1. ,92. ,18340. ,984000. ,. . . 

51200000. ,E( 1-5)=1. ,60. ,19510. ,558300. ,74368000. ,. . . 

F(l)=2750. ,G(1-3)=1. ,42. 2,2750. ,H(1-9)=1. ,40. 652,. . . 

2949. 22,1048754. 7,1. 86426E+08,3. 9622E+09,!. 17521E+11,. . . 

2. 12662E+12, 1. 10017E+14,L( 1-11)=1. ,53. 12,30506. 757,. . . 

1378867. , 1. 83946E+08,5. 2509E+09,3. 42441E+11,2. 85168E+12,. . . 

1. 43994E+14,!. 06498E+13 , -3. 42095E+14 ,K( 1-2)=13. 3,39. 9,. . . 

M(l-2) = 1. ,30. ,P(1)=25. ,Q(1-2) = 1. ,100. ,Y(1) = 1. ,Z(1) = 1. ,. . . 

S(l-2)=10. ,50. ,T(1-4)=1. ,4. ,8. ,0. 

PARAM R=l. 

*Changing zeta and omgn changes the specifications of your system 

CONST ZETA=0. 40, 0MGN=6. 0, 1C=0. 0 

INITIAL 

*This part of the program generates the standart time response 
T0TERR=0. 0 

C1=SQRT(1. 0-ZETA*ZETA) 

C2=0MGN*ZETA 

C3=0MGN*C1 

PHI=ATAN(C1/ZETA) 

* ASSIGN INITIAL GUESS FOR FUNCTION MINIMIZATION 
•'This part gives the starting point to the program 

IF (FLAG. GE. 0. ) GO TO 10 
A1=A10 

IF (FLAG. GE. 0. ) GO TO 10 
A2=A20 

IF (FLAG. GE. 0. ) GO TO 10 
B1=B10 

IF (FLAG. GE. 0. ) GO TO 10 
B2=B20 
10 CONTINUE 

*Put your variable coefficients here so program can change it 
Y(2)=A1 
Y(3)=A2 
Z(2)=B1 
Z(3)=B2 

*Write normal simulation program for your system 
DERIVATIVE 
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X0=TRNFR(8,10,0. ,H,L,X2*2. 849) 

X2=TRNFR( 0,2,0. ,F,G,X3) 

X3=TRNFR(4,4,0. ,D,E,X4) 

X4=TRNFR( 0,1,0. ,A,B,X6) 

X6=15.*ERG 

ERG=ER-X9 

X9=TRNFR(2,2,0. ,Y,Z,X0) 

ER=R-X0 

DYNAMICS 

*Standart time response is generated by this formula 
STD=1.0-(EXP(-C2*TIME)/C1)*SIN(C3*TIME+PHI) 

ERR=X0-STD 
ERRSQ=ERR*ERR 
T0TERR=INTGRL(0. ,ERRSQ) 

*As a cost function integral of error square is used 
TERMINAL 

*This part finishes the calculation when min value is reached. 

ERRFN = TOTERR 
CONTRL FINTIM=2. ,DELT=0. 01 
END 
STOP 
FORTRAN 

*From now on just declare your step values and amount of variables 
* MAIN PROGRAM FOR FUNCTION MINIMIZATION 
IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION X(4) ,STEP(4) ,0(4) ,QQ(4) ,W(4) 

C X(1)=K,X(2)=P INITIAL VALUES IN DSL ROUTINE 
STEP(1)=0. 5D0 
STEP(2)=5. ODO 
STEP(3)=0. 5D0 
STEP(4)=10. ODO 
N=4 

ITMAX = 100 
CFTOL = 1. OD-6 
ALPHA = 2. ODO 
BETA = 0. 5D0 
IPRINT = 0 
MINMAX = -1 

CALL H00KE( X , STEP , N , ITMAX , CFTOL , ALPHA , BETA , 

* CF,Q,QQ,W, IPRINT, MINMAX) 

STOP 

END 
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